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For crucial tests between theory and experiment, ab initio close coupling cal- 

O 

culations are carried out for photoionization of several oxgyen ions: O II, O 

£T) | III, O IV, O V. The relativistic fine structure and resonance effects are stud- 

O 

ied using the R-matrix and its relativistic variant the Breit Pauli R-matrix 
(BPRM) approximations. Extremely detailed comparison is made with high 

H , 

| resolution experimental measurements carried out in three different set-ups: 



Advanced Light Source at Berkeley, synchrotron radiation experiment with 



| a new modulator at University of Aarhiis, and synchrotron radiation experi- 



ment at University of Paris-Sud. The comparisons illustrate physical effects in 
photoionization such as (i) fine structure, (ii) resolution, and (iii) metastable 
components. Photoionization cross sections api of the ground and a few 
low lying excited states of these ions obtained in the experimental spectrum 
include combined features of these states. The general features of the mea- 
sured api were predicted in the earlier calculations in LS coupling. However, 
due to higher resolution achievement in the recent experiments, the theo- 
retically calculated resonances need to be resolved with much finer energy 
mesh for precise comparison. In addition, prominent resonant features are 
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observed in the measured spectra that can form from transitions allowed with 
relativistic fine structure, but not in LS coupling. Present results include 
calculations from both LS coupling and the relativistic BPRM method to de- 
cipher the features in photoionization spectra, a pi are obtained for ground 
and metastable (i) 2s 2 2p 3 ( 4 S° 2 D° 2 P°) states of O II, (ii) 2s 2 2p 2 ( 3 P, 1 D, 1 S) 
and 2s2p 3 ( 5 S°) states of O III, (iii) 2s 2 2p{ 2 P°j) and 2s2p 2 ( 4 Pj) levels of O 
IV, and (iv) 2s 2 ( 1 S') and 2s2p( 3 P°, 1 P°) states of O V. Most of the computed 
theoretical resonance structures are seen experimentally, although there are 
some discrepancies. One of the main conclusions of the present study is that 
resonances in photoionization cross sections of the ground and metastable 
states can be a diagnostic of experimental beam composition, with potential 

applications to astrophysical and laboratory plasma environments. 
PACS number (s): 32.80.Fb 
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I. INTRODUCTION 



Photoionization cross sections of many atoms and ions were studied in detail with reso- 
nance structures for ground and many excited states under the Opacity Project (OP, [1]). 
The close coupling (CC) R-matrix method, as employed under the OP, accurately considers 
the numerous autoionizing resonances in cross sections along interacting (overlapping) Ry- 
dberg series. The study under the OP resulted in new features in the cross sections, such 
as prominent PEC (photo-excitation of core) resonances. The radiative work under the OP 
has been extened to collisional and radiative processes under the Iron Project (IP, [2]) with 
inclusion of relativistic effects using the BPRM method. However, while the methods used 
are state-of-the-art, possible deficiencies are: (i) the resonance structures in the OP/IP cross 
sections are not always fully resolved, and (ii) the important near-threshold resonances may 
not be accurate owing to inadequate eigenfunction expansions. A prime need therefore is to 
benchmark the vast amount of theoretical OP/IP data that is now being utilized for global 
applications such as calculation of stellar opacities, photoionization rates in astrophysical 
models. 

On the experimental side, accurate photoionization cross sections are being measured 
at several places, such as the merged ion-photon beam facility at the University of Aarhiis 
(e.g. for C II [3]), the Electron Cyclotron Resonant Ion Source (ECR) at the University of 
Paris-Sud (e.g. for Xe V - Xe VIII [4]), and the Advanced Light Source (ALS) at Berekely 
(e.g. for O II [5]). Each of these three sophisticated experimental set-ups has different 
advantages in measurements, and the results often complement each other. Cross sections 
for O II have been measured at ALS for near threshold features in o P j of metastable states 
[5], including a wider photon energy range [6], and at the University of Aarhiis [7]. Detailed 
photoionization cross sections have recently been measured for the ground and low excited 
states of multiply charged oxygen ions for the first time for O III, O IV, and O V, at the 
synchrotron radiation source at University of Paris-Sud [8]. Study of oxygen ions is of great 
importance because of their abundances in terrestrial and planetary atmospheres and in 
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astrophysical objects. 

Comparison of all these results is extremely crucial to evaluate the accuracy of both the 
theoretical and experimental results, and to study the resonance phenomena. Theoretical 
work can also provide spectroscopic identification of resonances, and a quantitative mea- 
sure of the mixture of states in the experimental beam. The resonant structures and the 
background shapes of the recently measured photoionization cross sections of oxygen ions 
are found to be in close agreement overall with the some ealier theoretical cross sections [9] . 
However, certain features in the measured cross sections are not seen prominently or are 
missing in the calculated cross sections. 

One of the aims of this study is to study relativistic effects and fine structure. We 
distinguish these in terms of relativistic mixing and fine structure splitting with reference to 
pure LS coupling results. While relativistic effects per se might not be of great importance 
for these light ions, fine structure splitting manifests itself in many experimentally observed 
features. In addition, relativistic mixing via El transitions (e.g. singlet-triplet, doublet- 
quartet) results in additional resonances not allowed in LS coupling. The present study 
investigates the LS coupling and fine structure effects together with resonance phenonmena 
in detail. Thus to begin with it is necessary to do LS coupling calculations at high resolution. 
Comparison with experiment then reveals whether fine structure is also to be considered in 
order to analyze experimental data. In such cases the relativistic BPRM method is employed 
to compute cross sections. 

Detailed comparison of experimental and theoretical resonance positions and shapes also 
highlights the importance of correlation effects in theoretical calculations. Resonances may 
be shifted with respect to their respective thresholds of convergence, indicating missing 
correlation configurations in the wavefunction expansion. This again is crucial since the 
inherent accuracy of the close coupling calculations depends on the completeness of the 
eigenfunction expansions, in particular short-range correlation at low energies. 

Another motivation for these studies is that the R-matrix method, as developed for the 
OP/IP work, has been extended to a theoretically self-consistent treatment of photoion- 
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ization and unified electron-ion recombination including radiative and dielectronic recom- 
bination in an ab initio manner [10-12]. The total recombination cross sections have been 
benchmarked for accuracy by comparing with several experimental cross sections [11,12]. 
But while recombination cross sections require total contributions of photorecombination 
cross sections from all bound levels, obtained from the correnpoiding a PI using detailed 
balance, comparison of photoionization cross sections with experiment involves individual 
features of a P i of single states. For instance, photoionization cross sections including rel- 
ativists effects in the close coupling approximation, as developed under the IP, was first 
benchmarked by comparing with the measured a pi of C II [13]. 

We outline the theory and calculations in Sections II and III respectively, followed by 
a discussion of results for the four ions O II - O V in Section IV, and the conclusions in 
Section V. 

II. THEORY 

Present calculations are carried out in an ab initio manner using the close coupling R- 
matrix method (e.g. [14,15]) as in the Opacity Project [1] and the Iron Project [2]. The 
package of codes from the OP [15] have been extended to include relativistic effects in the 
Breit-Pauli approximation in the IP work [2,16,17] under the name Breit-Pauli R-matrix 
(BPRM) method. A brief discussion of the method is given below. 

In the close coupling approximation the total wavefunction of the (N+l)-electron ion is 
represented by the wavefunctions of the N-electron core, multiplied by the wavefunction of 
the outer electron: 

^ E (e + ion) = A^2xi(ion)9i + ^2c j § j (e + ion), (2.1) 

i j 

where Xi is the target (or the core) wavefunction in a specific state SiLiiii or level J(Ki and 
9i is the wavefunction for the (JV+l)-th electron in a channel labeled SiL^J^ikfl^ Jit); 
kf is its incident kinetic electron energy. $j's are the correlation functions of the (N+l)- 
electron system that account for short range correlation and the orthogonality between the 
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continuum and the bound orbitals. The complex resonant structures in photoionization are 
included through channel couplings. 

Relativistic effects are included through Breit-Pauli approximation in intermediate cou- 
pling. The Breit-Pauli Hamiltonian, as adopted in the IP work [2], has the following terms, 

n N+l - n N+1 + n N+1 + n N+l + n N+l , 

where H^+i is the nonrelativistic Hamiltonian, 

N + l ( cyy N+l 9 ) 

i=l { ' * j>i 1 V ) 

and the additional terms represent the one-body mass correction, the Darwin, and the spin- 
orbit interaction terms, respectively: 

* i * i 1 i i 'i 

The set of SLir are recoupled to obtain (e + ion) states with total Jtc, following the 
diagonalization of the (N+l)-electron Hamiltonian to solve 

H* P +1 V E = EV E . (2.5) 

Substitution of the wavefunction expansion introduces a set of coupled equations that are 
solved using the R-matrix approach. The continuun wavefunction, typ, describe the scatter- 
ing process with the free electron interacting with the target at positive energies (E > 0), 
while the solutions correspond to pure bound states $g at negative total energies (E < 0). 

The photoionization cross section (api) is proportional to the generalized line strength 
(S), 

opi = -^-—ujS. (2.6) 

where S is the generalized line strength: 

S L = | <^||D L ||^> | 2 , (2.7) 

D is the dipole operator; = J2i r i m "length" form, and Dy = — 2J2i m "velocity" 
form, where the sum corresponds to number of electrons. 
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III. CALCULATIONS 



Photoionization cross sections of O II, O III, O IV, and O V are obtained using wavefunc- 
tion expansions for residual core ions as developed in earlier works (Nahar [9]). Calculations 
are carried out both in non-relativistic LS coupling and in the BPRM approximation. 

The wavefunctions of the target or core ion were obtained from atomic structure cal- 
culations using a scaled Thomas-Fermi model in the code SUPERSTRUCTURE [18]. The 
computation of O II photoionization cross sections is repeated with the same wavefunction 
expansion of 12 terms of O III of configurations 2s 2 2p 2 , 2s2p 3 , 2s 2 2p3s, and 2p A as in Ref. 
[9]. The wavefunction expansion for O III consists of configurations 2s 2 p, 2s2p 2 , 2p 3 , 2s 2 3s, 
2s 2 3p, 2s 2 3d, 2s2p3s, and 2s2p3p of O IV. However, to include the fine structure in calculat- 
ing a pi for 2s2p 3 ( 5 S°) state of O III, an expansion of 15 levels of O IV with configurations 
2s 2 p, 2s2p 2 , and 2p 3 was used. The O IV eigenf unctions consist of a 20-level expansion of 
the core ion O V with configurations 2s 2 , 2s2p, 2p 2 , 2s3s, 2s3p, 2s3d. The wavefunction 
expansion for O V consists of 9 states of the core ion O VI with ls 2 2s, ls 2 2p, ls 2 3s, ls 2 3p, 
ls 2 3d, ls 2 4s, ls 2 Ap, ls 2 4d, and ls 2 4/. The calculated target level energies were replaced 
by the slightly different observed energies to obtain somewhat more accurate positions for 
resonances that belong to Rydberg series converging on to the target thresholds. 

The short-range bound channel correlation term in the CC expansion, the second term 
in Eq. (1), includes all possible (N+l)-electron configurations as in Ref. [9] for LS coupling 
calculations. BPRM calculations are carried out for O III and O IV to elicit fine structure 
details apparent in experimental data. For the BPRM calculations of O III, all configurations 
are considered with maximum orbital occupancies 2s 2 , 2p A , 3s 2 , 3p 2 , 3d 2 , while the K-shell 
remains filled. For O IV, the (N+l)-electron configurations are up to Is 2 , 2p 3 , 3s 2 , 3p 2 , 3d 2 , 
As 2 , Ap 2 , while K-shell is open with a single occupancy. 

Computations of erpj are carried out using the package of BPRM codes from the Iron 
Project [17]. The present calculations are considerably longer than earlier works for detailed 
comparison with experiments. They are repeated many times over, scanning with finer en- 
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ergy meshes as needed, to search for and delineate narrow resonances. For complex ions, 
resonances are extensive and overlapping, and are easy to miss when energy mesh is not fine 
enough. A direct comparison of theoretical and experimental resonant structures may not 
be possible if the resonances are narrower than the monochromatic bandwidth of the exper- 
imental beam. The resonances are usually convolved with a Gaussian distribution function 
with full width half maximum (FWHM) equal to the energy bandwidth of the monochrom- 
eter in the experiment. The average cross section, a™ I (E i ), at energy Ei convoluted over 
the FWHM, E Q , or the monochromatic bandwidth, is obtained as 



where x = E ~ Et . Present cross sections have been convolved with typical energy bandwidths 
used in various experiments for comparison. The photon energies are shifted slightly to 
match the measured ionization thresholds of each state exactly. 



Photoionization cross sections of the ground and several lowest excited states of O II, O 
III, O IV, and O V are obtained and compared with the recent precise and high-resolution 
measurements at the Advanced Light Source (ALS) at Berkeley, the University of Aarhiis, 
and the University of Paris-Sud. The resonances of detailed photoionization cross sections 
have been convolved with experimental bandwidth, as mentioned above, using a Gaussian 
distribution function to compare with the measured cross sections. Both the detailed and 
convolved opi are presented for each ion. 

Photoionization cross sections in LS coupling in general show very good agreement with 
experiment, indicating that relativistic effects are not significant for most states considered, 
especially for O II and O V. However, inclusion of relativistic effects is crucial for some 
metastable states for O III and O IV. 

The main objective of this report is to compare the theory and experiment for the exis- 
tence and the positions of resonant structures and their identification in the cross sections. 




(3.1) 



IV. RESULTS AND DISCUSSIONS 
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However, the exact shape and peak of the resonances, as well as the background cross 
sections, depend on (i) resolution, (ii) distribution function for the beam shape, and (hi) 
fractions or mixture of states in the beam. The third point is difficult to determine exactly 
as the beam composition may vary with energy during measurements. With some numeri- 
cal variations, the theoretical and experimental shapes can be brought in closer agreement, 
which however is not the main emphasis of the present work. The Gaussian distribution 
function, which peaks around the center, is used for convolution of theoretical cross sec- 
tions. This may not describe exactly the distribution function of the experimental beam. 
Experimental calibration of the residual ion threshold energies is also an important factor 
for locating resonance positions. 

Observed resonant features in photoionization of oxygen ions have been identified in Refs. 
[5,7,6,8]. Resonances can be identified using quantum defect analysis: E r = E t —z 2 /u 2 , where 
E r is the resonance energy, E t is the next core threshold energy in Rydbergs, z = Z — n c 
is the ion charge of the target where n c is the number of core electrons, and v = n — /i is 
the effective quantum number with \i being the quantum defect. \i is large for low angular 
momenta, such as for an s-electron, and decreases with higher angular momenta approaching 
near zero for d orbitals in oxygen ions. 

The BPRM results have many more resonances than in LS coupling due to (i) an incre- 
ment in number of Rydberg series of resonances with increased number of core thresholds 
with fine structure and (ii) relativistic mixing of El allowed and intercombination transitions. 
LS coupling results are approximately the statistical average over the fine structure. 

Results for each ion are presented and discussed in separate subsections below. 

A. Photoionization of O II 

Photoionization cross sections of the ground 4 S° state and the metastable 2 P°, 2 D° states 
of O II in the low energy region, 30 - 36 eV, are shown in Fig. 1. The figure has seven panels: 
the lowest three panels, Fig. la,b,c, present the detailed a PI of 4 S , ° (red curve), 2 D° (cyan), 
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2po (g reen ) states respectively. These cross sections include extensive resonances due to 
Rydberg series of autoionizing states. They become extremely narrow as the autoinizing 
states converge on to the core thresholds, 1 D ) 1 P, 5 S°, etc. and their integrated strength 
is usually small. They are wiped out during convolution, and are not observed, as can be 
seen in the convolved cross sections of 4 S°, 2 D°, and 2 P° states in panel d. The colors 
of curves identify the resonances of the states they belong to. Panel e presents sums the 
convolved cross sections with mixing fractions 15% for 2 P°, 4% for 2 D°, and 43% for 4 S°, 
as deduced from experimental beam mixture in the ALS experiment. The summed cross 
sections can be compared directly with the a P j measured at ALS (panel b), and at Aarhiis 
(panel a). Although the mixing proportions are somewhat different in the two experiments 
the prominent features appear in both. The ALS cross sections have much higher resolution 
than those from Aarhiis because of the much narrower bandwidth. However, the ALS 
values are relative while the Aarhiis cross sections are absolute. Very good agreement can 
be seen between the theoretical and experimental features. The most significant point is 
that all the resonant features in the observed photoionization spectra have now been found 
in the calculated cross sections. Previous calculations [9] used much coarser energy mesh in 
calculating a PI and missed some resonances while others were not adequately resolved. 

Fig. 2 presents a P j of the same states, the ground 4 S , ° and the metastable 2 P°, 2 D°, 
but in the high energy region, 36 - 46 eV. It also has seven panels: the lowest three panels 
(a,b,c) present the detailed a PI of 4 S° (red curve), 2 D° (cyan), 2 P° (green) states. The 
extensive narrow resonance structures of highly excited Rydberg states converging on to the 
various core thresholds are nearly wiped out during convolution, as seen in the convolved 
cross sections for 4 S°, 2 D°, and 2 P° states together in panel d. Identical colors are used 
for both the detailed and the convolved cross sections of each state. Panel c presents sum 
of the convolved a PI with mixtures of 15% for 2 P°, 4% for 2 D°, and 43% for 4 S° as in 
the ALS experimental beam. The summed cross sections are compared with the measured 
a PI at ALS (panel b) and Aarhiis (panel a). Very good agreement can be seen for most of 
features theoretically and experimentally. While most of the observed resonant features in 
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the photoionization spectra are seen in the calculated cross sections, some discrepancies can 
be noted. The observed resonance at 37.5 eV, belonging to 2 D° state, is shifted to the right 
(higher energy) by about 0.5 eV in the calculated spectrum. This may have resulted from 
some missing correlation configurations. Two resonances around 37.6 and 38.9 eV are also 
not found in the calculated cross sections. As these resonances were not seen in the high 
resolution ALS data, it is possible that they belong to some other states in the mixture. 

The resonance series have been identified in Refs. [5,7,6]. For example, the resonances 
below 31.8 eV (threshold for 2 D°) in Fig. 1 are from 2p 3 { 2 P°) state which can photoionize 
to 2 S, 2 P, and 2 D states with 2p electron going to ed or es continua. Hence, the resonances 
belong to 2p 2 ( 1 D)nd( 2 S 2 P 2 D) and 2p 2 ( 1 D)ns( 2 D) Rydberg series of threshold 2p 2 ( 1 D). 
The resonances of higher angular momentum in general dominate over the lower ones. The 
nd series can be differentiated from ns series via quantum defects and similarity in features. 
For example, of the two distinct resonance series seen in the low energy region of Fig. Id, the 
first series (with higher peaks) actually represents three overlapping nd series of resonances, 
2p 2 ( 1 D)nd( 2 S 2 P 2 D) and the second one (lower peaks) is the ns series, 2p 2 ( 1 D)ns( 2 D). 
The table of identified resonances are given in Refs. [5,7]. 

B. Photoionization of O III 

The first measurement of photoionization cross sections for O III is reported by the Uni- 
versity of Paris-Sud [8]. To compare with the experiment, we repeat calculations to obtain 
high-resolution cross sections for the three lowest states, 2s 2 2p 2 ( 3 P, 1 D, 1 S) of the ground 
configuration, and the metastable 2s2p 3 ( 5 S°) state of O III using the 23CC wavefunction 
expansion for the residual Boron-like ion O IV of Nahar [9]. The detailed cross sections of 
2s 2 2p 2 ( 3 P, 1 D, 1 S) states are presented in three different panels, a-c, of Fig. 3 ( 3 P - red, 1 D 
- blue, and 1 S - green) and the convolved cross sections are in panel d and the weighted 
sum in panel e. Present detailed cross sections with higher resolution display much more 
narrow resonances compared to the earlier calculations. Comparison between theory and 



11 



experiment, panels e and f, shows that all the resonances observed for the three states are 
accounted for in the calculated values. 

While relativistic effects are not important to bring out the main features of the three 
2s 2 2p 2 ( 3 P, 1 D, 1 S) states, a closer examination of the experimental results shows that there 
are additional features, especially in the low energy region (see panel f). We therefore 
consider the first odd parity metastable state 2s2p 3 ( 5 S°) of O III and include fine structure 
through BPRM calculations. Prominent resonant features are found in the low energy region 
of <jpj of this state that are not allowed in LS coupling, but are due to relativistic channel 
mixing. In particular the very first resonance at ~48 eV in the experimentally observed 
cross sections is due to photoionization of the high-lying metastable 2s2p 3 ( 5 S°). 

We outline the differences in a PI between the LS coupling and the BPRM fine structure 
cross sections in Fig. 4. Figs. 4a and b present (Tpj( 5 S°) in LS coupling and in the BPRM 
approximation respectively. The cross sections in the two panels are very similar from a 
photon energy of 56.4 eV onwards. Whereas a PI is zero below this threshold in LS coupling, 
inclusion of relativistic effects in the BPRM calculations shows the region filled with extensive 
resonances with almost zero background cross section. At 56.4 eV, 2s2p 3 ( 5 S°) of O III ionizes 
leaving the core ion in the excited 4 P state of O IV, which is the first threshold to the 
electron continuum. In LS coupling, 2s2p 3 ( 5 S°) state photoionizes only into a 5 P (e + ion) 
angular and spin symmetry. The 5 P symmetry can not be formed from the ground state 
continuum of O IV i.e. 2 P°nl. In other words, 56.4 eV is the ionization threshold for the 5 S° 
metastable state of O III, and equals the sum of 47.4 eV, for ionization into the ground 2 P° 
state, and 8.85 eV, for excitation of the core to the 2s2p 2 ( 4 P) threshold. No features can 
be formed below this threshold in LS coupling. However, fine structure allows ionization of 
the 2s2p 3 ( 5 5 < 2) level of O III into the core ground state of O IV through relativistic mixing. 

Expanded and detailed a PI of O III below the 4 P threshold of O IV are presented in 
Fig. 4c to enable comparison with experiment. Patterns of various resonance complexes, 
as denoted by R™ series in the figure, can be discerned. The resonances belong to several 
Rydberg series of autoionizing levels, 4 Pi/2,3/2,5/2nd( 5 Pj) and 4 Pi/2,3/2,5/2^s( 5 Pj). The 5 Pi, 2 ,3 
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levels are allowed to decay to 2 P°l/2, 3/2np( 1,3 Pi j2 ,3) levels via radiationless transitions. 
Decay to other levels, such as 1,3 -D li2) 3, 1,3 i 7 2,3 are also possible through some relativistic 
mixing. However, these transitions are very weak. 

The convolved cross sections with experimental monochromatic bandwidth for 5 S° state 
are shown Fig. 3d. Many theoeretically computed resonances are difficult to observe in 
the experimental data shown in the top panel of Fig. 3 since their width is less than the 
monochromatic bandwidth of the detector in the experiment. A weighted sum of the ground 
and all metastable states included is given in Fig 3e. The excellent agreement between 
theory and experiment shows that the theoretical O III cross sections accurately reproduce 
the corresponding features in the experimental data. 

C. Photoionization cross sections of O IV 

Fine structure is considered in calculating photoionization cross sections of O III to 
show that large number of resonances can be delineated. With increasing ion charge, and 
depending on the ionic species, fine structure can manifest itself more prominently which is 
the case for photoionization of boron-like O IV. We present opi for fine structure levels of 
the ground LS term, 2s 2 2p( 2 P°/ 2 3/ , 2 ), and of the first metastable term, 2s2p 2 ( 4 P 1 / 2 ,3/2,5/2) of 
O IV. Earlier works on C II, an isoelectronic partner of O IV, showed that some resonances 
observed in the experiment were missing from the LS coupling results [3]. Both relativistic 
mixing and fine structure splitting were found to play a role in photoionization of B-like ions 
and explained the missing resonances [13]. 

Fig. 5 presents relativistic apj for the levels 2 Pjy 2 3/2 ari< ^ 4 -^1/2,3/2,5/2 • Cross sections for 
the fine structure levels of each term show very similar features. The number of resonances 
in the BPRM results is much larger in comparison to those in LS coupling, as expected 
owing to fine structure splitting of target thresholds resulting in a larger number of Rydberg 
series of resonances. For example, the 12 LS terms of the core ion O V correspond to 20 fine 
structure levels, which provide 8 more series of resonances. Some addtional resonances are 
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also produced due to to relativistic fine structure couplings (A J = 0, ± 1) not allowed via 
LS coupling. 

Figs. 5c-e show prominent resonant features, although with almost zero background, in 
the low energy BPRM cross sections of 4 P levels. Extensive narrow resonances are seen to 
converge at about 78.7 eV, confirming enhancement in the observed photoionization cross 
section (Fig. 6c). Similar to the 5 S° state of O III, channel couplings of fine structure 
levels are important in photoionization of the metastable 4 P term of O IV in forming these 
resonances. Without fine structure mixing there will be no cross section below 78.7 eV, the 
ionization threshold in LS coupling for 4 P state leaving the core O V in 3 P° state. At this 
LS term energy, both the BPRM and LS coupling calculations give the same photoionzation 
cross sections (analagous to the 5 S° state of O III). However, fine structure levels coupling 
lowers the ionization threshold from 78.7 eV by about 10.2 eV to 68.5 eV; where 10.2 eV is 
the energy difference between the target thresholds 2s 2 ( 1 S') and 2s2p( 3 P°) of O V (the O 
IV 4 P couples and photoionizes into the latter state). These low energy resonances can be 
explained as follows. The 2s2p 2 4 P state can photoionize to 4 D°, 4 P°, 4 S° SLn symmetries of 
the (e + ion) continua. Interaction of channels with fine structure allows couplings with total 
angular momenta (SL)Jvr Jit = 1/2° and 3/2°. The autoionizing levels 2s2p( 3 P°)nd( 4 D°) 
thereby decay to Jn=l/2° and 3/2° continua 2s 2 ( 1 S , )ep( 2 P°) via radiation-less transitions. 

To compare these photoionization cross sections directly with the experimentally mea- 
sured ones, presented in Fig. 6a, the cross sections are convolved using a Gaussian distribu- 
tion function with FWHM equal to the monochromatic beamwidth. The opi of each level 
are then statistically averaged (£i (2Si + 1)(2L; + l)a PI /J2i (2$ + 1)(2L, + 1)) to obtain 
the convolved cross sections of 2 P° and 4 P states of O IV. These convolved and averaged 
a pi are presented in Fig. 6c (red - 2 P° and green 4 P). 0>/( 4 P) are multiplied by 9% and 
o PI P° by 91% before addition to obtain the weighted sum a PI . These weighted cross sections 
are given in Fig. 6b. The top panel presents the experimentally measured photoionization 
cross secitons at the University of Paris-Sud [8] . 

Resonances in the energy region below the ground 2 P° ionization threshold, 77.4 eV (Fig. 
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6), are due to the metastable 4 P state (green curve in Fig. 6c). The metastable state is 
relatively easy to excite as it corresponds to a A n = transition (2s 2 2p — 2s2p 2 ) and its 
features are seen. The first resonance in Fig. 6 is at about 69.7 eV and its energy position 
agrees both theoretically and experimentally. However, some shift in positions (by ~ 0.5 eV) 
can be seen for the next three resonances. The exact reason for this shift is not clear, but 
a possible reason could be the effect of correlations of configuration interactions for the 4 P 
symmetry. For the rest of the resonances the measured positions agree with the theoretical 
ones. One exception is the measured resonance at about 79.5 eV which is not found either 
in the 2 P° or the 4 P calculated cross sections. 

The calculated peaks of resonances in Fig. 6 are seen mostly higher than the measured 
ones. These peaks may be reduced via radiative decays, or dielectronic recombination, if 
the radiative rates are comparable to autoionization rates. However, the radiative decay 
rates of the relevant excited core thresholds 10 8 , 10 9 sec" 1 ) of O V are several orders of 
magnitude lower than typical autoionization rate (~ 10 13-14 sec" 1 ), and thus should cause 
very little damping, i.e., reduction in peak values. As mentioned earlier, the calculated 
peaks are affected by resolution or different energy meshes, Gaussian energy width, and 
percentage abundances of states in the beam. Particular mixture fractions can even eliminate 
resonances, as seen in the comparison between Fig. 6 (b) and (c). Number of resonances in 
convolved cross sections in Fig. 6c are vanished in the weighted sum in Fig. 6b. Hence, for 
identification of resonances, the experimental spectrum should be compared with the pure 
convolved cross sections, i.e. between Fig. 6a and c, and as identified in [8]. 

The most significant feature in the measured cross sections is the prominent resonance 
at about 77 eV (Fig. 6a). This is a combined effect from the near threshold resonances of 
the 2 P° ground state and converging resonances of 4 P levels to core threshold 3 P°/ 2 , as seen 
in Fig. 5. Although narrow resonances are not often observed becasue of low strengths, 
the correspondence in energy positions indicates that the metastable 4 P photoionization is 
making a contribution to the feature. The difference between the theoretical and experimen- 
tal strengths in the feature is an indicative that the mixtures of states may not be exactly 
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the same as quoted in the experiment. However, the fact that Fig. 6b exhibits the features 
seen experimentally is of potential importance in the determination of fractional abundances 
of the ground and excited states in the ion source. The background cross section can be 
reduced or enhanced in differenct energy regions by taking various fractions. The fractions 
91% and 9% used in the present simulation may not exact match the experimental ones. 
Also, theoretical resonances are resolved with a much finer energy mesh than the experi- 
ment beam width, and hence appear sharper. The Gaussian convolution does not broaden 
the features as in the experiment, which indicates that the experimental beam distribution 
function was not purely Gaussian. In addition, the FWHM may not remain constant at all 
energies; a convolved peak can be flattened with a larger bandwidth. 

Comparison between the theoretical and experimental photoionization cross sections in 
Fig. 6 highlights two of the main points in this work: (i) a mixture of the ground 2 P° and 
metastable 4 P states in the experimental beam represents the observed features; both theory 
and experiment show the same resonances, and (ii) the detailed cross sections depends on 
the fractional mixture of ionic states. 

D. Photoionization cross sections of O V 

Photoionization cross sections of O V are obtained for the ground 2s 2 ( 1 1 S) and the excited 
2s2p( 3 P°, 1 P°) states. The detailed features in opi are shown in Fig. 7a-c: 3 P° - red, 1 D - 
green, 1 P° - blue. The resonances are more resolved than those reported earlier [9]. 

The cross sections of the states are convolved with experimental bandwidth to compare 
with measured values at the University of Paris-Sud [8]. The convolved cross sections are 
presented in Fig. 7d, keeping the original colors for each state for easy identification. For 
comparison with experiment, thes convolved cross sections are summed with weighting fac- 
tors 75% for 1 S, 23% for 3 P°, and 2% for 1 P°, whereas the fractions in the experimental 
beam as reported in Ref. [8] are 65% for the ground state, 1 S, and 35% for excited 3 P° state. 
The weighted sum in Fig. 7e is compared with the measured spectra in Fig. 7f. The reason 
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for inclusion of 2s2p( 1 P°) cross sections is similarity of features. Although the measured 
features in the cross sections below 103.6 eV, the ionization threshold for 2sp2( 3 P°) state, 
appear as small kinks in Fig. 7f, they correspond to the resonant features appearing at the 
same energy positions of 2s2p( 1 P°) state of O V. 

A clear rise seen in the background cross sections at about 103.6 eV corresponds to the 
3 P° ionization threshold, and at about 113.85 eV to 1 S ionizaton threshold, as noted in 
7e and 7f. However, differences in the background cross sections are found. Since present 
mixture includes higher percentage of the 1 S ground state, similar enhancement is expected 
in the weighted average. The calculated resonant peaks are also higher than the observed 
ones. Discussion regarding the differences in peaks is made above. However, the main 
objective is met, that is, all observed resonant features, shapes and energy positions are 
reproduced in the calculated spectra. Identifications of the resonances are given in Ref. [8]. 

V. CONCLUSION 

Photoionization cross section of oxygen in several stages of ionization have been obtained 
for the ground and low lying excited states. These cross sections are compared with measured 
cross sections at three different experimental set-ups. 

While the experimental measurement includes combined features of resonances belonging 
to different Rydberg series of states, and of a mixture of different states, theoretical work 
can evaluate them for individaul states. 

Certain features in the measured cross sections are not seen prominently or are missing 
in previous calculations. Present report shows that the missing features are due to (i) lower 
resolution of previous calculations or (ii) due to relativistic fine structure effects, not included 
in the previous calculations. For most cases of this relatively low Z element, LS coupling 
calculations are often adequate. However, for some metastable states, fine structure is crucial 
for the low energy resonances which allowed only the couplings of fine structure levels. It is 
the fine structure components, rather than the relativistic effects, that is important for the 
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formation of these high peak resonances with nearly zero backgound. Inclusion of relativistic 
effects will certainly be more exact and accurate way to carry out the calculations. 
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FIGURES 
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FIG. 1. Photoionization cross sections apj of ground 4 S" 5 and metastable 2 P°, 2 D° states of 
II in low energy region, 30 - 36 eV: Panels a, b, c - detailed a pi of 2 P° (green), 2 D° (cyan), 
and ^S (red) states respectively; d - the same cross sections, but are convolved with experimental 
monochrormeter bandwidth: e - sum of the convolved cross sections weighted by experimental 
beam mixtures to compare with measured opi\ f , g - measured opi by Covington et al. (2001) and 
by Kjeldsen et al. (2002), respectively. 
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FIG. 2. Photoionization cross sections apj of ground 4 S° and metastable 2 P°, 2 D° states of 
II in the higher energy region. 36 - 46 eV: Panels a, b, c - detailed a pi of 2 P° (green), 2 D° (cyan), 
and 4 5" 5 (red) states respectively; d - the same cross sections, but are convolved with experimental 
monochrormeter bandwidth: e - sum of the convolved cross sections weighted by experimental 
beam mixtures to compare with measured apr, f, g - measured opi by Aguilar et al. (2003) and 
by Kjeldsen et al. (2002), respectively. 
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FIG. 3. Photoionization cross sections a pi of III: Panels a, b, c - 2s 2 2p 2 ( s P, 1 D, 1 S) states; 
d - convolved cross sections of 3 P, 1 D, 1 S, 5 S°; e - weighted sum of convolved cross sections; f - 
experimental cross sections of the combined states at University of Paris-Sud by Champeaux et al. 
(2003). 
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FIG. 4. Photoionization cross sections a pi of excited metastable 2s2p 3 ( 5 S*2) state of III: 
Panel a - in LS coupling; b - in BPRM approximation; c - 2s2p 2 A Pnl{R) series identification! in 
low energy region. 
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FIG. 5. Photoionization cross sections api of IV in BPRM approximation: a, b - P° 
levels of ground configuration, 2s 2 2p; c, d , f - excited 2s2p 2 ( 4 P 1 / 2 ,3/2,5/2) levels. 
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FIG. 6. Comparison of photoionization cross sections apj of IV with experiment: c - red 
and green curves are statistical sum of convolved cross sections of levels 2 P°/ 2 3/2 an d 4 -Pi/2.3/2,5/2 
respectively; b - weighted sum of the cross sections; c - measured cross sections at University of 
Paris-Sud by Champeaux et al. (2003). 
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FIG. 7. Photoionization cross sections apj of V: a, b, c - Detailed opi of ground 2s 2 ( 1 S') 
and excited 2s2p( 3 P°, 1 P°) states; d - convolved cross sections, e - weighted sum of convolved cross 
sections; f - measured cross sections at University of Paris-Sud by Champeaux et al. (2003). 
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